Demixing can occur in binary hard-sphere mixtures with negative non-additivity 
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A binary fluid mixture of non-additive hard spheres characterized by a size ratio 7 = C2/C1 < 1 
and a non-additivity parameter A = 2oi2 / (01 + 02) — 1 is considered in infinitely many dimensions. 
From the equation of state in the second virial approximation (which is exact in the limit d — > 00) a 
demixing transition with a critical consolute point at a packing fraction scaling as n ~ d2~ d is found, 
even for slightly negative non-additivity, if A > — | (hi7) 2 . Arguments concerning the stability of 
the demixing with respect to freezing are provided. 
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The crystallization of a hard-sphere fluid, first ob- 
served in computer simulation []] and at that time con- 
troversial, is a clear and presently well established exam- 
ple of an entropy-driven phase transition. Yet, up until 
now it has defied a rigorous statistical mechanical proof. 
Other phase transitions governed by entropy, in partic- 
ular the phase separation in binary hard-core mixtures, 
has only been proven [2( for a simple two-dimensional lat- 
tice model of two types of particles. In the absence of an 
exactly solvable model in three dimensions, which could 
provide insight into the detailed mechanisms leading to 
phase separation in athermal systems such as hard-core 
mixtures, different strategies have been adopted to ad- 
dress this problem. For instance, to consider a geometry 
that leads to a non trivial volume-driven phase separa- 
tion, as in the case of a mixture of parallel hard cubes 
0. Or else^to explore the situation in higher spatial 
dimensions 



Due to the interest (both theoretical and in prac- 
tical applications) of phase separation, the demixing 
problem in additive hard-sphere mixtures has received 
a lot of attention in the literature. An analysis of the 
solution of the Percus-Yevick equation for binary ad- 
ditive hard-sphere mixtures [f| led to the conclusion 
that no phase separation into two fluid phases existed 
in these systems. The same conclusion is reached if 
one considers the most popular equation of state pro- 
posed for such mixtures, namely the Boublfk-Mansoori- 
Carnahan-Starling-Leland (BMCSL) [6j equation of 
state. For a long time the belief was that this was a true 
physical feature. Nevertheless, this belief started to be 
seriously questioned after Biben and Hansen obtained 
fluid-fluid segregation in additive hard-sphere mixtures 
out of the solution of the Ornstein-Zernike equation with 
the Rogers- Young closure and subsequent work has con- 
centrated on attempting to clarify the issue. Coussaert 
and Baus Q have proposed an equation of state with im- 



* andres@unex.es Jhttp^/ww^unex^es / eweb / f i s t e or / and res/] 
1 malopez@servidor.unam.mx http: //miquiztli . cie .unam.mx/xml/tc 
On sabbatical leave from Centro de Investigacion en Energia 
U.N.A.M., Temixco, Morelos 62580 (Mexico) 



proved virial behavior for a binary additive hard-sphere 
mixture that predicts a fluid-fluid transition at very high 
pressures (metastable with respect to a fluid-solid one). 
On the other hand, Regnaut et al. have examined 
the connection between empirical expressions for the con- 
tact values of the pair distribution functions and the ex- 
istence of fluid-fluid separation in mixtures of additive 
hard spheres. Further, in the case of highly asymmetric 
binary additive hard-sphere mixtures, the depletion ef- 
fect has been invoked as the ph ysic al mechanism behind 
demixing (see for instance Ref. U£j and the bibliography 
therein). Finally, demixing in mixtures of additive hard 
spheres has been examined recently j 1 1 1 using the low 
density expansion of the pressure by adding successively 
one more exact virial coefficient (up to the sixth virial 
coefficient). In this latter work it was found that al- 
ready within the second virial coefficient approximation 
the fluid separates into two phases of different composi- 
tion with a lower consolute critical point. 

In contrast to the above results, which have the draw- 
back of having been derived under various approxi- 
mations and are therefore open to question and con- 
troversy, non-additive hard-core systems with positive 
non-additivity are certainly known to exhibit fluid-fluid 
demixing, although again this has not been rigorously 
proved in general. The celebrated Widom-Rowlinson 
model represents a prototype system that allows the 
detailed study of such a phase transition, an aspect that 
continues to be of interest in the recent literature [l3[ . 
Provided fluid-fluid segregation really occurs in additive 
hard-sphere mixtures, where size asymmetry would be 
the source of the transition, it is not unreasonable to 
expect that, given a certain degree of (high) size asym- 
metry, demixing may also be present in the case of hard- 
sphere mixtures with small negative non-additivity. This 
feature, however, seems to have hardly received any at- 
tention [l4j |. The purpose of this Communication is to 
address this problem and provide evidence in favor of 
the statement posed in the title of the paper. To do so, 
we will not work in three-dimensional space, but rather 
consider the limit of infinitely many dimensions in which 
nnr rpsijlt will be exact. 

/ft/mlh/ 1 — 1 i—j 

Although there had been a few earlier papers [lj, 
dealing with hard spheres in dimensions greater than 
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three, it was after the pioneer work of Frisch et al. 
[l7| in which they showed that the classical hard-sphere 
fluid in infinitely many dimensions was amenable to 
full analytical solution, that studies of high dimensional 
hard-sphere systems became common over the years 

he ii! nisi mm. Thc fact that features such as 

the freezing transition are present in all dimensionalities 
(except for d = 1) and the parallel between high spa- 
tial dimensions and limiting high density situations that 
seems to exist in fluids (with the added bonus of greater 
mathematical simplicity as one increases the number of 
dimensions) suggest that one can gain insight into the 
thermodynamic behavior of say three-dimensional sys- 
tems by looking at a similar problem in higher dimen- 
sions. As a matter of fact, the very elegant work of 
Carmesin et al. |2lj has exploited this approach to il- 
lustrate phase separation in a hard-sphere mixture with 
positive non-additivity in infinite spatial dimensionality. 
We now consider the more general case also for d — ► oo 
in which the non-additivity may take negative values. 

Let us consider a binary mixture of non-additive hard 
spheres of diameters o\ and a 2 in d dimensions. The hard 
core of the interaction between a sphere of species 1 and a 
sphere of species 2 is o\ 2 = + ct 2 )(1 + A), where the 
parameter A characterizes the degree of non-additivity 
of the interactions. Further assume (something that will 
become exact in the limit d — > 00 [2l| ) that the equation 
of state of the mixture is described by the second virial 
coefficient only, namely 



p = pk B T [1 + B 2 {xi)p] 



(1) 



where p is the pressure, p is the number density, k B is 
Boltzmann's constant, T is the temperature, and 



B 2 (xi) = v d 2 d ~ 1 (x\ai + x\a d + 2x 1 x 2 a d 2 



(2) 



is the second virial coefficient, x\ and x 2 = 1 — x\ being 
the mole fractions and v d = (7r/4) d / 2 /r(l + d/2) being 
the volume of a (i-dimensional sphere of unit diameter. 
The Gibbs free energy per particle is (in units of fcgT) 



where B 2 = B 2 /2 d - 1 v d af, A; = (Ai/ax) d , and A d = 
d/2 d ~ 1 Vd- Next we take the limit d — > 00 and assume 
that the volume ratio 7 = 7 d is kept fixed and that 
there is a (slight) non-additivity A = d~ 2 A such that 
the scaled non-additivity parameter A is also kept fixed 
in this limit. Thus, the second virial coefficient can be 



approximated by B 2 = B, 



(0) 



5(0) _ /_ 
B 2 - [Xx 



Y 1 ' 2 ) and B?> = 



B^ 

(i) 



0(d 2 ), where 
z 1 / 2 ^, with 



K = i (In 7) + 2 A. Let us remark that in order to find 
a consolute critical point, it is essential to keep the term 
of order <i _1 if A < 0. The equation of state (J5J can 
then be inverted to yield y = j/ ' + y^d~ x + 0(d~ 2 ), 

with y<® = JpB® and y« = -\ (l - y^B^/B^' 



In turn, the Gibbs free energy © becomes g — g^d + 
g^ + Oid- 1 ), with = 2i/(°), g^ = ^ Xi In (x l \ l ) + 

In (^Ady^ / B^^j + 2y( 1 \ while the chemical potentials 

9 - xi {dg/dx 1 ) p 

h 0(d~ 1 ), where 



Mi = 9 + x 2 {dg/dxi) p and p 2 



are given by pn 



N d 



Mi 



(o) 



2p 1/2 , Mi X) =hx[A dXl Xi yJp/B, 



i/V4 0) 



{x 2 /x 1 ){jp) 1 / 2 B { 2 1) /B 2 :o) and M2 is obtained from Ml by 
the changes x\ <-> x 2 , \\ — > A2/7, 7 — » I/7, p — > pj, 
B 2 -» fla/7- 

The coordinates of the critical point are readily found 
to be 



Xlc = 



7 



3/4 



1 + ^3/4 



Pc = 



(1+7 1/4 ) Z 
47^ 2 



(7) 



Note that x\ c is independent of A. The coexistence 
curve, which has to be obtained numerically, follows from 
the conditions p,f~ \xa,P) — P^\xb,p) (i = 1,2) where 
x\ = xa and X\ — xb are the mole fractions of the coex- 
isting phases. Once the critical consolute point has been 
identified in the pressure/concentration plane, we can ob- 
tain the critical density. The dominant behavior of B 2 



g = x\ In (x\pkf) + x 2 In (x 2 pA 2 ) + 2B 2 (xi)p, (3) at the critical point is B 2 {x\ c ) — 7/ (l 



7 1/4 + 7 



l/2\ 



where Aj (i = 1, 2) are the thermal de Broglie wave- 
lengths. Given a size ratio 7 = <j 2 jo\ < f, a value 
of A, and a dimensionality d, the consolute critical point 
(xi c ,p c ) is the solution to (d 2 g/dx 2 ) p = (d 3 g/dxf) p = 0, 
provided of course it exists. Then, from Eq. one can 
get the critical density p c . 

We now introduce the scaled quantities 

p = 2 d - 1 v d d- 2 paf/k B T, y = <r x B 2 p. (4) 

Consequently, Eqs. and © can be rewritten as 

p = y(y + d- 1 )/B 2 , (5) 



Xi In {x l \ l ) + In (A d y/B 2 ^ + 2dy, (6) 



while y, 



(0) 



(f + 7 1 / 4 ) /2 (l — 7 X /4 + ^1/2) Hence, 



the critical density readily follows after substitution in 
the scaling relation given in Eq. For our purposes 
it is also convenient to consider the packing fraction 
defined as r\ = Vdpuf {x\ + x 2 ^f) and its scaled version 
rj = d~ 1 2 d rj j23|. At the critical point, this latter takes 
the nice expression 



Vc 



(7 



1/8. 



•7" 



-1/8V 



(8) 



Figure n shows x\ Cl p c , and rj c as functions of 7 and 
in the two latter cases for A = —0.1 (negative non- 
additivity), A = (additive mixture) and A = 0.1 
(positive non-additivity). The previous results clearly 
indicate that a demixing transition is possible not only 
for additive or positively non-additive mixtures but even 
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FIG. 1: (Color online) Plot of X\ c (inset, upper panel), p c 
(upper panel), and rj c (lower panel) as functions of 7. The 
(scaled) critical pressure and packing fraction are displayed 
for three different values of the (scaled) non-additivity pa- 
rameter: A = —0.1 (solid lines), A = (dash-dotted lines), 
and A = 0.1 (dotted lines). Note the non-monotonic depen- 
dence of p c and rj c on 7 (for the latter case, see inset in the 
lower panel). 

for negative non-additivities. The only requirement is 
K > 0, i.e. A > — |(ln7) 2 or, equivalently, A > 
— i(ln7) 2 . The curve representing the threshold situ- 
ation A = — I (ln7) 2 is plotted in Fig. [3 where we have 
also displayed A, as obtained from Eq. @, as a function 
of 7 for three different values of the critical packing frac- 
tion: rj c = 1, rj c = 1.5, and rj c = 2. These choices for 
rjc are meant to be illustrative and have been taken after 
the following considerations. 

One may reasonably wonder whether the demixing 
we have obtained for negative non-additivity will occur 
for packing fractions within the stable fluid regime and 
where the equation of state is well represented by the sec- 
ond virial approximation. A natural way to look into this 
issue would be to compare with the close-packing value 
rj cp . Unfortunately, rj cp is not known in the case of mix- 
tures. Nevertheless, some insight about it can be gained 
by examining the parallel case of a one-component fluid 
in infinitely many dimensions. In such a case, there 
exist known upper bounds for rj cp |22|. Further, another 
(lower) estimate can be given by taking the contribu- 
tions of the second and third virial coefficients to be of 
a similar order of magnitude. While this is of course not 
conclusive, all these estimates for rj cp may be shown to 
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FIG. 2: (Color online) Plot of A vs 7 for different values 
of the critical packing fraction: rj c = 2 (thin solid line), 
rjc = 1.5 (dash-dotted line), and rj c = 1 (dotted line). The 
lowest thick solid line corresponds to the threshold condi- 
tion A = — I (hi7) 2 . A demixing transition with the scaling 
properties of the text is only possible for mixtures with non- 
additivities and size ratios represented by points above the 
thick curve. On the other hand, a transition with a (scaled) 
critical packing fraction smaller than a given value of rjc is 
only possible for points above the corresponding curve. The 
threshold curve goes to minus infinity as 7 — > but the other 
curves have a minimum that depends on the choice of rj c (see 
inset). This implies that for extremely asymmetric mixtures 
(7 —* 0) demixing at a given finite value of rj c becomes possi- 
ble only if the non-additivity is positive and sufficiently large. 



diverge as d — > 00, suggesting that the fluid- fluid phase- 
separation may indeed take place. Provided the (scaled) 
packing fraction at freezing rjf is different from zero, the 
demixing transition may be stable and not preempted 
by a fluid-solid transition. Again rjf is unknown but we 
may once more recur to the one-component case. For 
this system, Colot and Baus have conjectured that 
(rjf /rj cv ) 1 / d becomes independent of d for high d. Fur- 
ther, from the analysis of the results in d — 3, 4, 5 [l6Ll20| . 
and d = 7 |2^| one finds that rjf w 1.3. Since at freez- 
ing or melting the Hclmholtz free energies of the fluid 
and the solid should be of the same order of magnitude, 
by considering the former given by the second virial ap- 
proximation and the latter as obtained from free volume 
theory with the estimate (?7//?7cp) 1/d ~ 0.8, we obtain the 
rough estimate rjf m 2.3. Irrespective of the numbers, the 
point is that these results seem to confirm that rjf is dif- 
ferent from zero and finite. Therefore, even if the range 
of values of negative A in which stable demixing occurs 
is limited and restricted to highly asymmetric mixtures, 
as indicated in Fig. |2 the important issue is that it is 
certainly there. So the question of whether demixing can 
occur in binary mixtures of hard spheres with negative 
non-additivity can be given a positive answer. 

While the high dimensionality limit has allowed us 
to address the problem in a mathematically simple and 
clear-cut way, the possibility of demixing with negative 



4 



non-additivity is not an artifact of that limit. Demix- 
ing is known to occur for positive non-additive binary 
mixtures of hard spheres in three dimensions and com- 
pelling evidence in the additive case exists, at least in 
the metastable fluid region. Even though in a three- 
dimensional mixture the equation of state is certainly 
more complicated than Eq. and the demixing transi- 
tion reported here for negative non-additivity is possibly 
metastable with respect to the freezing transition, the 
main effects at work (namely the competition between 
depletion due to size asymmetry and hetero-coordination 
due to negative non-additivity) are also present. In fact, 
it is interesting to point out that Roth et al. |14|. using 
the approximation of an effective one-component fluid 
with pair interactions to describe a binary mixture of 
non-additive hard spheres and employing an empirical 



rule based on the effective second virial coefficient, have 
also suggested that demixing is possible for small nega- 
tive non-additivity and high size asymmetry. Our exact 
results lend support to this suggestion based on approx- 
imate calculations and confirm the fact that the limit 
d — ► oo allows one to get a caricature or toy model to 
highlight features already present in real systems. 
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